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Performance of quantum process estimation is naturally limited to fundamental, random, and systematic 
imperfections in preparations and measurements. These imperfections may lead to considerable errors in the 
process reconstruction due to the fact that standard data analysis techniques presume ideal devices. Here, by 
utilizing generic auxiliary quantum or classical correlations, we provide a framework for estimation of quan- 
tum dynamics via a single measurement apparatus. By construction, this approach can be applied to quantum 
tomography schemes with calibrated faulty state generators and analyzers. Specifically, we present a general- 
ization of "Direct Characterization of Quantum Dynamics" [M. Mohseni and D. A. Lidar, Phys. Rev. Lett. 97, 
170501 (2006)] with an imperfect Bell-state analyzer. We demonstrate that, for several physically relevant noisy 
preparations and measurements, only classical correlations and small data processing overhead are sufficient to 
accomplish the full system identification. Furthermore, we provide the optimal input states for which the error 
amplification due to inversion on the measurement data is minimal. 



PACS numbers: 03.65.Wj, 03.67.-a, 03.67.Pp 
I. INTRODUCTION 

Quantum measurement theory imposes fundamental limita- 
tions on the information extractable from a quantum system. 
Although the evolution of quantum systems can be described 
deterministically, the measurement operation always leads to 
nondeterministic outcomes. In order to obtain a desired ac- 
curacy, measurement of a particular observable needs to be 
repeated over an ensemble of identical quantum systems. In 
addition, for systems with many degrees of freedom, one usu- 
ally needs to measure a set of non-commuting observables 
corresponding to independent parameters of the system. Char- 
acterization of state or dynamics of a quantum system can 
be achieved by a family of methods known as quantum to- 
mography llIlQl - In particular, quantum process tomography 
provides a general experimental procedure for estimating dy- 
namics of a system which has an unknown interaction with its 
embeddiiig environment for discrete or continuous degrees of 
freedom Q HI B HI 01 ■ In these methods, the full information 
is obtained by a complete set of experimental settings associ- 
ated with the set of required input states and non-commuting 
measurements. In recent developments HI 0, S HI S it 
has been demonstrated that the minimal number of required 
experimental settings can indeed be substantially reduced by 
using degrees of freedom of auxiliary quantum systems corre- 
lated with the system of interest. 

Generally it is possible to completely characterize a quan- 
tum device with a single experimental setting. A correlated 
input state of the combined system and an ancilla is subjected 
to the unknown process, and a generalized measurement or 
Positive Operator Valued Measure (POVM) is performed at 
the output 1 1 1, 12J. However, in order to realize such a gener- 
alized measurement one needs to effectively generate many- 
body interactions lfT2ll which are not naturally available and/or 



controllable. Quantum simulation of such many-body interac- 
tions is in principle possible, but generally requires an expo- 
nentially large number of single- and two-body interactions 
with respect to system's degrees of freedom. An alternative 
method that circumvents the requirement for many-body in- 
teractions, yet allows simultaneous non-commuting observ- 
ables through a single measurement setting, is known as Di- 
rect Characterization of Quantum Dynamics (DCQD) HSl- 
The construction of the full information about the dynami- 
cal process is then possible via preparation of a set of mu- 
tually unbiased entangled input states over a subspace of the 
total Hilbert space of the principle system and an ancilla HI ■ 
The DCQD approach was originally developed with the as- 
sumptions of ideal (i.e., error-free) quantum state preparation, 
measurement, and ancilla channels. However, in a realis- 
tic estimation process, due to decoherence, limited prepara- 
tion/measurement accuracies, or other imperfections, certain 
errors may occur hindering the overall procedure. 

In this work, we introduce an experimental procedure for 
using generic two-body interactions to perform quantum pro- 
cess estimation on a subsystem of interest. We employ this ap- 
proach to generalize the DCQD scheme to the cases in which 
the preparations and measurements are realized with known 
systematic faulty operations. We demonstrate that in some 
specific, but physically motivated, noise models, such as the 
generalized depolarizing channels, only classical correlation 
between the system and ancilla suffice. Moreover, for these 
situations the data processing overhead is fairly small in com- 
parison to the ideal DCQD. Given a noise model, one can 
find the optimal input states by minimizing the errors incurred 
through the inversion of experimental data. Thus, we provide 
the optimal input states for reducing the inversion errors in the 
noiseless DCQD scheme. 

The structure of the paper is as follows. In Sec. [Ill we 
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set the framework for process tomography schemes where we 
have faulty — rather than ideal — faulty Bell-state analyz- 
ers, emphasizing the DCQD approach. Next, in Sec. HUl we 
demonstrate the applicability of our framework through some 
simple, yet important examples of noise models. We conclude 
by summarizing the paper in Sec.lIVI 



II. CHARACTERIZATION OF QUANTUM PROCESSES 
WITH A FAULTY BELL-STATE ANALYZER 

Let us consider a given quantum system composed of two 
correlated physical subsystems A and B. For a time duration 
At the two subsystems are decoupled, thus experiencing dif- 
ferent quantum processes, and then they interact with each 
other again. The task is to characterize the unknown quan- 
tum process acting on the subsystem of interest A, assuming 
we have prior knowledge about the dynamics of subsystem B 
plus their initial and final correlations. Another similar sce- 
nario can also be envisioned. Given two controllable quantum 
systems A and B that are made to sufficiently interact before 
and after a time duration At, we wish to estimate the unknown 
dynamics acting on system A for such time interval, assum- 
ing the dynamics of the ancilla system B and the interaction 
between two systems is known with a certain accuracy. 

Much progress has been made in creating and characteriz- 
ing two-body correlations in a variety of physical systems and 
interactions, including nuclear magnetic resonance (NMR) 
systems interacting through an Ising Hamiltonian together 
with refocusing or dynamical decoupling techniques ifisll . 
atoms/molecules in cavity quantum electrodynamics (QED) 
ifTill . trapped ions interacting via the Jaynes-Cummings 
Hamiltonian driven by laser pulses and vibrational degrees 
of freedom ITsll . and photons correlated in one or many de- 
grees of freedom, e.g., generated by parametric-down conver- 
sion |[l6ll or four-wave mixi ng f lTll. Other approaches include 
spin-coupled quantum dots lIlSll . superconducting qubits lfl9ll 
controlled by external electric and/or magnetic fields, and 
chromophoric complexes coupled through Forster/Dexter in- 
teractions and monitored or controlled via ultra-fast nonlin- 
ear spectroscopy ll20ll . However, in almost all of these sys- 
tems, the entangled Bell-state preparations and measurements, 
which generically are the basic building blocks of quantum 
information processing, hardly achieve high fidelities; at the 
very least they are certain to be imperfect at some level, and 
this will limit their use for tomography. Our goal is to deter- 
mine the optimal states and measurement strategy that will 
minimize the deleterious effects of the nonidealities — as- 
sumed known — on process tomography. 

We consider the cases where we can simulate initial or fi- 
nal two-body correlations in the above schemes by performing 
an ideal (generalized) Bell-state preparation (or measurement) 
followed by a known faulty completely -positive (CP) quantum 
map acting on the entire systems involved. It should be noted 
that not every CP maps can be written as concatenation of two 
other CP maps. In other words, there exist CP maps that are 
"indivisible", in the sense that, for such a map T, there do not 
exist CP maps Ti and T2 such that T — T2T1, where neither 
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noisy BSP noisy BSM 

FIG. 1: (color online) Schematic of a faulty DCQD, with imper- 
fect or noisy Bell-state preparation (BSP) and Bell-state measure- 
ment (BSM). 



Ti or T2 are unitary ll2lll . Nonetheless, all full -rank CP maps 
— in the sense of the Kraus representation — are divisible. 

Here, we include quantum maps acting on system B in the 
preparation or measurement maps. This approach naturally 
provides a generalization of the DCQD scheme to the cases of 
faulty preparations, measurement and ancilla channels where 
the noise is already known — see Fig. [T] For simplicity, in 
this work we concentrate only on one-qubit systems and the 
DCQD scheme (summarized in Table H]). However, general- 
ization of the framework is straightforward for other process 
estimation schemes and for DCQD on qudit systems with d 
being a power of a prime (according to the construction of 
Ref. i]). 

Let us consider the qubit of interest A and the ancil- 
lary qubit B prepared in the maximally entangled state 
\^+)ab = (|00) + \n))AB/V2. We first apply a 
known quantum error map f^'^ to A and B: (p) = 
Epqr. Xpgrscr^f po-f erf , where p = |$+)($+| and 
{(Jo = 1, o"!, (72, fa} are the identity and Pauli operator for 
a single qubit. Next, we apply the unknown quantum map 
£ to qubit A; this is what we are trying to determine: 

Finally, we apply a known quantum error map before the Bell- 
state measurement. Note that in this approach any error on 
the ancilla channel can be absorbed into either £^^^ or f ^^'^ 
The total map acting on the combined system AB is then 



TABLE I: The ideal direct characterization of single-qubit x- Here 
l$i> = a|00) + (3\n), \^t)..(y) = «l + +).iy) -).(^) 
where |q| / |/?| / and Im(dP) / 0, and {|0), |1)}, {|±)4, 
{|±)j;} are eigenstates of the Pauli operators cr^, a^, and ay. 
is the projector on the Bell state \^^), and similarly for the other 
projectors. See Refs. Determination of optimal values of a 

and (3 is discussed in the text. 



input state 


Bell-state measurement 


output Xmn 






Xoo, Xii> X22, X33 


*i> 




X03 I X12 




±p^+,p^- ±p^- 


XOl, X23 




p^+ ±p^_,p^_ ±p^,+ 


X02 , X13 
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£(T) =£(f) o£o£«, given by 



yh'j ~ / ^ XmnXpqrsXp'q'r's' ^ 



E 



where the parameters Xp\rs and xl'vi-'s' known (from 



(f) 



p' q' r s 



calibration of the operational/systematic errors in the prepa- 
ration and measurement devices). By defining oJmnp's' = 

(_l)(5,„o-l)(5,.„-l)5,„,. + (5,.o-l)(5,,o-l)5„., p™„ = 



l^pp'qq'rr'ss' ^rnnp' s' Xpqrs Xp' q'r' s'"p' " q' "p "q P"r "s "r'"s 

we have 



A — A 
Xmn'^mPmn'^n ■ 



By construction, the parameters Xi'lrs and Xp'q'r's' 
priori known, as are the matrices p,„„ which are functions of 
and the initial state p. Therefore, in order to develop 
a generalized DCQD scheme for the systems with faulty Bell- 
state preparation (BSP) and measurement (BSM), we need to 
do it for a set of modified (input) states rather than a pure 
Bell-state type input. Expanding in the Bell basis yields 



where the A{x^'\ X^^K {p^'"*}) matrix contains full informa- 
tion about all faulty experimental conditions. Given x^'^ 
x'"^\ and the standard DCQD input set {p*^*^}, one can cal- 
culate the A matrix. The standard DCQD experimental data 
(analysis) will also determine Ix*-^-*). Now, if the A matrix 
is invertible, from Eq. ([T]i one can obtain x by inversion: 
\x) = A^^lx^^)). The invertibility of the A matrix, namely 
dct A 7^ 0, depends on the input states {jO*-*'} and the noise 
operations x^'^ and x'*^). It may happen that the A matrix 
becomes ill-conditioned ll22ll for a specific set of input states 
(for some given noise operations x*^'^ and x*^^'')- In such cases, 
even small errors (whether operational, stochastic, or round- 
off) in estimation of x'^'' can be amplified dramatically after 
multiplication by A^^. This in turn may render the estimation 
of X (the sought-for unknown map £) completely unreliable. 
To minimize the statistical errors, the input states should be 
chosen such that dct A is as far from zero as possible. There- 
fore, the optimal input states {p^pt} [optimal in the sense of 
minimizing statistical errors] for given x''-* and x^'^'' are ob- 
tained via maximizing det A. A similaT faithfulness measure 
has already been used in Refs. 10, 123I1 . In Appendix iBl we 
derive the optimal input states for the case of the ideal DCQD 
scheme. 



E\fefc 
kk' ^mn^ 



kk' iztkk 



where A: 



kk' 



Tr[P''*'p„„], P 



kk' 



>0 



B^) (B 



>k' 



and 



\B'') for k = 0,1,2,3 corresponds to the Bell-states |$+), 
Ivp"), and |$"), respectively, where |<I>='=) = (|00) ± 
|ll))/\/2, = (|01) ± \10))/V2. (Henceforth through- 

out this manuscript, superscripts refer to the Bell-state basis 
and subscripts refer to the Pauli operator basis.) The Af^Jj^s are 
known functions of uimnp's', x''^ x'''^ and p. Therefore, the 
overall output state can be rewritten as follows: 



III. PROCESS ESTIMATION WITH SPECIFIC NOISY 
DEVICES 



In the following we describe several examples which de- 
scribe relevant physical noise models. 



A 



kk' 



kk'mn ''^mn^^^'-' m 



We now apply the standard DCQD data analysis to esti- 
mate the matrix elements of x'^-* (representing £^'^'>). After 
performing a BSM, i.e., measuring {P-'-'}^j^Q on this state, we 
obtain the Bell-state \B^) with probability 



kk'mn ^^kk'mn^''^'^'' 



(3) 

kk'. 



X'^„Tt[P" a^P'''' a^]. Although this ex- 



where A 

pression can be made more compact by using Pauli identities, 
the current form is convenient for our purposes. 

A similar set of equations for the standard DCQD inputs 
{p^^^i=o can also be written. We represent all of these equa- 
tions in a compact vector form as 



X<^^)-A|X), 



(1) 



A. Depolarizing channels: Correlated noise 



An important and practically relevant example is the situ- 
ation in which 5^'' and S'^^^ both are two-qubit (hence corre- 
lated) depolarizing channels pt^l 



where e and e' could be independent of each other or corre- 
lated (e.g., £ — e'). These errors result in the following noisy 
data processing of the measurement results of DCQD: 



Tx[£{p^'^)P" 



{l-e){l-e') 
16 



Tr[£-(]l)(g) 1] 



-Tr[£(l) ® IP"\ 



e(l-e') 



Tr[£(p('))] + ee'TT[£{p^'^)P"]. (2) 
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For the Hamiltonian identification task ll26l IztIi . E{p) = 
^-iHt p^iHt (vvhichis unital: E{t) ~ t, and trace-preserving: 
Tr[£(p)] = 1), we obtain 

Tr[f ^ £e'Tr[f + (1 - ee')/4. (3) 

This relation provides a simple connection between the ideal 
and the noisy data processing rules. Another feature of Eq. (|3]l 
is that it is valid irrespective of the values of e and e' (7^ 0). 
This implies that, whether e and e' are in the range which 
makes the noisy preparation/BSM separable or not ll28i I29I1 . 
the simplicity and applicability of (the modified) DCQD re- 
main intact. In other words, entanglement is not an imperative 
in the DCQD algorithm. 

A generalization of this noise model is the case in which the 
preparations are modified based on a generalized two-qubit 
depolarizing channels ll24ll : 




-1 (g) 1 +e[/p(*'C/t, 



in which U is an already known two-qubit unitary operator. 
To simplify the following discussion we assume that BSMs 
are noiseless = I). Finding the exphcit form of x'*' is 
straightforward. We use the form 

where denotes the constrained summation in which the 
case (m, n) — (0, 0) has been excluded. Using the identity 

OkOi(Tk = (-l)^"'''"cr;, we have: \ Yll.b=o^o. ® ^bP^^a ® 
CTh = /0+fl®l,or equivalently: Ylb^QPabi^a ® crbpcra » 
(Jb = 1 (g 1, where pab = 1/3 except for pqq = —1. In 
addition, we expand U in the {am <E) (Jn}mn=o basis: U = 
J2mn O'mnO-m <8) cr„. Altogether, these relations yield 

(p) = '~4~ X! PmnCr,n ® cr„pcr„ (g) (T„ 

mn 

mn.m'n' 

Hence, we obtain: Xmnmn = Pmn(l - e)/4 + £|a,„„p (the 
diagonal elements) and xtLm'n' = for (m, n) ^ 

(to', n') (the off-diagonal elements). In a compact form, the 
effect of this noise channel can be expressed as follows: 

Tr[f (pW)FJJ] ^ 

^i-^^fc^Tr[f (1) ® 1] + £^il_£lTr[f (1) ® IP"] 

+ ^^^-j^Tv[£iUp^'^U^)] + ee'Ti[£iUp^'^U^)P"]. (4) 

Under trace-preserving and unitality conditions, the final data 
processing is thus modified as follows: 

Tr[£:(p(*))pj'J'] ^ ee'Ti[£{Up'''^U^)P"] + (1 - ee')/4. (5) 



B. Depolarizing channels: Uncorrelated noise 

We assume that the input states and our measurements are 
diluted by depolarizing channels ll28il29ll acting separately on 
the principal and ancilla qubits, i.e., D (g)!), where V acts on 
a general single-qubit state p as follows: (p) = 1 + ep, 

or equivalently: = Yl^j=QPj'^jP'^j^ where po = (1 + 

3£)/4 and pi = P2 = = (1 ^ ^)/4, and positivity and 
complete-positivity of I?e require — l/3^e^l Isoll . 

As a special case we specialize on the characterization of 
the diagonal elements Xkk- This is particularly important in 
Hamiltonian identification tasks |26,,27J. It can be easily seen 
that for Bell-states P^^ we obtain 

pkkV^. 1^1® l+£2pfcfc. 

4 

Thus, to estimate Xkk, the necessary data processing is mod- 
ified as in Eqs. ^ and (|3]l by replacing: ee' {^s'Y and 
i (recall that p'"' = |<I>+)($+|). Here we have as- 
sumed that the input (measurement) depolarizing parameter 
is e {e'). This result implies that to estimate the diagonal el- 
ements Xkk, whether under correlated noise or uncorrelated 
noise, the DCQD scheme is robust and classical data process- 
ing is modified in a simple fashion. This has immediate appli- 
cations to the task of Hamiltonian identification ll26ll . 

C. Generalized depolarizing channels 

Here, we assume that the input states and/or measurements 
are diluted such that they effectively lead to (known) Bell- 
diagonal input states and/or Bell-diagonal measurements. 
Thus we obtain 

(i)8ev-3 

P ^ A^i'=0 tji'F 7 

■ B , t -I -I 

P33 _£> V f' P-J J 

This noise results in the following noisy data processing of the 
measurement results of DCQD: 

Tr[£(p(»))P""] ^ Y.vr £«:'4.Tr[£:(p(''))PJ'^"']. (6) 

That is, every measurement result of the new setting is a linear 
combination of the ideal results. If we define the vector |p) = 
{pij)^ , where p.y = Tr[£(p(*^)P^^], namely 

|p) ^ (Tr[£(p(0))p00],Tr[£(pW)P"],...,Tr[f(p(3))p33])^, 

and the matrix Aij^i'j' = eaie'^y, then (|6]) can be written as 
the following linear matrix transformation (see AppendixlAli: 

|p)^A|p). (7) 

If we arrange the output elements as in Table U we will have 

|p) = C|p), (8) 



Although this is not as simple as Eq. (|3]l, it yet retains a con- where C is the (constant) coefficient matrix, hence |p) 
siderable simplicity. AC " ^ | p ) . 
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IV. SUMMARY 

We have provided a scheme for utilizing auxiliary quantum 
correlations to perform process estimation tasks with faulty 
quantum operations. We have demonstrated our approach 
via generalizing the ideal scheme of Direct Characterization 
of Quantum Dynamics (DCQD) where the required prepara- 
tions and measurements could be noisy. It has been shown 
that when the systematic faulty operations are of the form of 
depolarizing channels, the overhead data processing is fairly 
simple. Moreover, these examples have revealed that for the 
DCQD scheme, entanglement is secondary. This, in turn, 
broadens the range of applicability of our scheme to quan- 



tum systems with certain controllable classical correlations 
of their subsystems. Therefore, our proposed method may 
have near-term applications to a variety of realistic quantum 
systems/devices with the current state-of-technology, such as 
trapped ions, liquid-state NMR, optical lattices, and entangled 
pairs of photons. 
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Mathemaical Sciences (PIMS), and the USC Center for Quan- 
tum Information Science and Technology for funding. 



APPENDIX A: EXPLICIT FORM OF EQ. © 



Table U suggests that if, instead of the conventional BSMs, we consider the expression 



Tr[P™£(p<°))] 
Tr[P"£(p(°>)] 

Tr[P='3£(p<">)] 



Tr[pO°£(p<i): 
Tr[P"£(p<i>; 
Tr[P™£(p 



(!)■ 



Tr[pii£(p(i' 
Tr[p0°f(p(2) 
Tr[P^^£(p(^) 
Tr[P™£(p(2' 
Tr[p33f(p(2) 
Tr[P™£(p<3> 
Tr[P"£(p<3) 
Tr[P™e(p(3) 
^Tr[P^^£(p(^' 



+Tr[p33£(p<i))] 
+Tr[P^=£(p(i))] 
-Tr[p33e(p(i))] 
-T¥[p22£(p(i))] 
+Tr[pii£(p(=^))] 
+Tr[P=*3f(p<^')] 
-Tr[P"£(p(=')] 
-Tr[p2=£(p(^))] 
+Ti'[P^^£(p<-^))] 
+Tr[p33£(p<3))] 
-Tr[p22e(p(3))] 
-Tr[piie(p(^')]/ 



1 1 

1 1 
1 -1 

1 -1 



1 1 
1 -1 



1 1 
1 1 

-1 1 / 



/Tr[P™£(p<°))] \ 
' Tr[P"£(pW)] \ 
Tr[P""£(p<0))] 
Tr[P^^£(p<''>)] 
Tr[P™f(p<i>)] 
Tr[P"£(p<i))] 
Tr[P^^£(p<i))] 
Tr[p33i;(p(i))] 
Tr[P™f(p(2))] 
Tr[P"f(p<^')] 
Tr[P""f(p<2))] 
Tr[P-'*-'*f (p<^>)] 
Tr[P™£(p(-'*>)] 
Tr[P"f(p<-''))] 
Tr[p22£(p<-''))] 

\Tr[P=»=»£-(p(-''))] 



the results of the measurements can be related to the Xmn elements in a more straightforward fashion. Here, the constant 
coefficient matrix (C) relates the results of the new BSMs (p) to those of the conventional BSMs (p). 



APPENDIX B: OPTIMAL INPUT STATES FOR THE IDEAL (NOISELESS) DCQD 



Here, we find the optimal input states for the ideal DCQD. The idea is to choose the input states such that the (linear) inversion 
on the experimental data (to read out x matrix elements) can be performed reliably. That is, the goal should be to make the 
coefficient matrix as far from singular matrices as possible. Maximizing the determinant of this matrix is a sufficient condition 
to guarantee its reliable invertibility, hence in turn, minimal error propagation. 

The data obtained from the measurements (BSMs) are Tr [£{p^'^^ )P^^ , where {p^*^} corresponds to the first column of Table|I] 
respectively, for i = 0, 1, 2, 3. We parameterize the input states by a = cos 6 and /3 = e**^ sin {ip ^ kn, k E Z). Using relation 
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FIG. 2: The value of | det A| vs 9 and ip. Here, the coefficient matrix A{8, tp) relates the experimental outcomes to the unknown elements 
of the superoperator Ix'^') = ■'^Ix)- The input states for the standard DCQD, as defined in Table|l] e.g., 1$+) = q|00) + are 
parameterized as a = cos 6 and /3 = e"^ sin 9 (ip ^ kn, k G Z). The optimal input states are associated to those parameters for which 
I det A I has its maximal value 1, leading to a minimal statistical error. 



(lAlb . one can express Eq. ([T]i, for the standard DCQD as the following: 



Tr| 
Tr| 
Tr| 

Tr[P™f (p< 
Tr[P"f(p< 
Tr[P°°S{p^ 
Tr[P"e(p< 
Tr[P™f(p< 
Tr[p22£(p( 
Tr[P™e(p' 
Tr[P=''''£(p' 
Tr[P™£(p< 
Tr[P"£(p' 
Tr[P™e(p' 
\Tr[P='^£(p' 



■P™£(p<i))] 
■P"£(p(i))] 

:p22f(p<i))] 
:p3^£(p<i>)] 

^')]+Tr[P=^^ 
^>)]+Tr[P^^ 
^')]-Tr[P^=* 
^')]-Ti-[p22 
^>)]+Tr[pi^ 
^')]+Tr[P^3 
^')]-Tr[pii 
^')]-Tr[P^^ 
''>)]+Tr[p22 
*>)]+Tr[P'*^ 
''')]-Tr[p22 
■*')]-Tr[P" 



e(p<=')] 
e(p<=')] 
£(P<^')] 
£(P<^')] 
£(p<")] 
e(p<=')] 
£(P<^')] 
e(p'^')] 

£(pW)] 

£(P<*')] 
e(p<*')] , 
£(P<*')]/ 
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where in the coefficient matrix A{9, (p) we have x 
matrix is obtained as 



cos 29, y = sin 29 sin (p, and z = sin 26 cos ip. The determinant of this 



IdetAI 



sin^ 49 sin^ ip, 



which attains its maximum value 1 at (6* = tt/8 + kn/i, p = tt/2 + k'n), Vfc, k' E Z — Fig. [2] Therefore, the optimal input 
states {p''^} for the standard DCQD are as in TableUin which fi and are either of the pairs calculated from the above maximal 
set of 6 and (p. A simple calculation shows that the amount of entanglement (exactly speaking, concurrence \3U ) of the optimal 
non-maximally entangled input states is 1/ \/2 (independent of (p). 
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